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We employ the density functional theory to study a structural transition of MnO from Bl (rock- 
salt) to B8 (NiAs) structures that was observed experimentally at pressures around 100 GPa. We 
utilize all-electron description as well as norm-conserving pseudopotentials and demonstrate that 
these two approaches can significantly differ in quantitative predictions. We explicitly show that 
even small-core pseudopotentials exhibit transferability inaccuracies for quantities sensitive to the 
energy differences between high- and low-spin polarizations of valence electrons. 
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I. INTRODUCTION 

Transition metal oxides have received a considerable 
attention of condensed matter physicists in the past sev- 
eral decades. One of the reasons for such an elevated 
interest is that the standard band theory methods, such 
as the density functional theory (DFT) represented by 
the local density approximation (LDA) or by its semi- 
local generalization (GGA), do not satisfactorily describe 
these materials. For example, the LDA/GGA predicts 
incorrect equilibrium crystal structure for FeOi^ In the 
case of equilibrium properties of MnO, these approxima- 
tions perform reasonably well as was shown by several 
groups. The source of deficiencies of simple band theo- 
ries is agreed to be underestimation of correlations among 
3c? electrons. Improved account for these correlations is 
needed to adequately capture Mott insulating mechanism 
and the competing charge transfer effects. 

Several methods have been devised to resolve the issue 
of strongly correlated electrons, such as GW, LDA+U, 
self-interaction correction (SIC) or hybrid exchange- 
correlation functionals that combine local exchange from 
LDA or GGA with non-local Fock exchange. Perfor- 
mance of several of these approaches has recently been 
thoroughly compared in Ref. [H, using the pressure- 
induced high-spin to low-spin transition 6 in MnO as a 
benchmark test. It turned out that although the cur- 
rent computational techniques work well for some as- 
pect of the electronic structure (such as the equilibrium 
volume or the character of the ground state at ambient 
pressure), they are not quite dependable if one considers 
high-pressure properties of this material. 

One of the goals of this paper is to provide additional 
benchmarks for the electronic structure of manganese 
oxide at high pressure, expanding thus the test set of 
Ref. H. Our main message does not directly relate to the 
3d electrons but to a rather surprising finding that the 
deep core states and the type of method used for their 
description can significantly influence certain properties 
generally considered to be related to the valence elec- 
trons only. An example of such a quantity is the critical 
pressure P c for the pressure- induced transition from Bl 



(rocksalt) to B8 (NiAs) structure that has been reported 
recently This transition is accompanied by a collapse 
of the local magnetic moments of Mn atoms& as well as 
with electronic transition to a metallic stated 

Another key motivation for this study of core-valence 
interactions in crystalline MnO is our interest in appli- 
cations of quantum Monte Carlo (QMC) techniques to 
transition metal oxides and other materials containing 
strongly correlated electrons. The accurate representa- 
tion of physical ions by norm-conserving pseudopoten- 
tials is essential for obtaining meaningful and predictive 
quantum Monte Carlo results and it is therefore crucial 
to assess the accuracy and limits of such pseudopoten- 
tial Hamiltonians. Ideally, one would compare results 
of pseudopotential QMC with all-electron QMC. Unfor- 
tunately, such a route is prohibitively computationally 
expensive at present. We have therefore decided to carry 
out a careful study within the DFT, Hartree-Fock (HF) 
approximation and combination of these methods, since 
the accuracy of one-particle orbitals is important in QMC 
as mentioned below and explained elsewherei 10 ! 11 Similar 
problems were inv estigat ed in several density-functional 
studies (e.g., Refs. fT3fl3ll 1 4il 1 5h and provided important 
signals that the choice of the methods and treatment of 
pseudopotentials can affect the final results. 

Our study provides new insights in several important 
and complementary directions. First, we employ Dirac- 
Fock pseudopotentials, which nominally should not re- 
quire any additional corrections (such as the DFT non- 
linear core correction as commented later on). The rea- 
son for this choice is that these pseudopotentials were 
designed to be accurate for correlated wave function cal- 
culations and, in fact, several tests have been carried 
out within both post-Hartree-Fock (configuration inter- 
action) and quantum Monte Carlo methods with quite 
satisfactory resultsi 10 ' 11 Second, we wanted to quantify 
the pseudopotential vs. all-electron differences within the 
DFT, HF and hybrid functionals since our previous^ 
(and current preliminary) calculations show that the hy- 
brid functionals provide the most optimal one-particle 
orbitals, which minimize the fixed- node error in the quan- 
tum Monte Carlo method. The use of two indepen- 
dent codes and basis sets provides a solid computational 
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framework for accurate benchmarking of these issues. In 
addition, study of high-pressure phase adds another im- 
portant test, since the dramatic changes in electronic 
structure probe the pseudopotentials in the regime, for 
which their reliability cannot be taken for granted. 



II. AMBIENT PRESSURE 

Before we proceed with investigation of the Bl — ► B8 
phase transition we first concentrate on energy differ- 
ences between various magnetic phases at ambient pres- 
sure. We apply the all-electron method side by side with 
the pseudopotential approach. For this purpose we have 
modified the WIEN2k package^ so that one can use 
either the original all-electron core solver or the norm- 
conserving pseudopotentials (PP). The valence states are 
expanded in the LAPW basis (extended by local orbitals) 
in both cases. This allows for a systematic comparison 
of full cores with PPs, since the only point where such 
calculations differ is the treatment of the core electrons. 
Similar combination of LAPW and pseudopotentials was 
already used to look into some aspects of PP applications 
to transition metals*^ 

In order to minimize errors associated with partition- 
ing core-valence correlations in pseudopotential calcula- 
tions, we have chosen pseudopotentials with 3s and 3p 
states in the valence space, so that the eliminated core 
corresponded to the Ne atom. We tested two types 
of norm-conserving PPs constructed from Dirac-Fock 
atomic solutions. The first one, denoted as YN, is a 
soft pseudopotential generated using Troullier-Martins 
construction^ The second PP, labeled as STU, is the 
so-called energy- adjusted pseudopotential, I 8 which has 
the —Z e ff/r attractive ionic part retained. 

In addition to performance of pseudopotentials we were 
interested in the impact of exact exchange on electronic 
structure of manganese oxide. To this end we have em- 
ployed the CRYSTAL2003 code±2 that uses Gaussian 
type atomic orbitals as a basis for one-particle wave func- 
tions. In the following we refer to this method as LCAO. 
The exact exchange was studied by means of two hybrid 
exchange-correlation functionals. We utilized B3LYP as 
well as somewhat simpler PBEO functional. The later 
one is schematically written as 
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where E PBE and E PBE are exchange and correlation 
parts of the PBE-GGA^ E BF is the exact (Fock) ex- 
change and a > 0. This type of hybrid mixing has the 
advantage that it directly relates to the PBE-GGA ap- 
proximation employed in our LAPW calculations. We 
use a rather small value a = 1/10 for the Fock mixing 
weight to account for the (nearly) metallic behavior at 
pressures we are going to investigate later. In the rest of 
the paper we abbreviate this functional as PBEOio- With 
the standard choice of a — 1/4 (Ref . l2lh we ran into se- 
rious difficulties in converging the Kohn-Sham equations 
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TABLE I: Energy difference A = Eafm - E NM per MnO at 
lattice constant 4.43 A together with local moments m formed 
on Mn atoms in the AFM phase. The disagreement between 
A aii and A pp has opposite signs in PBE-GGA and Hartee- 
Fock (HF) calculations and is considerably reduced for hybrid 
functionals, especially for B3LYP. 



when the lattice was compressed to the extent needed to 
reach the Bl — > B8 transition. Another alternative and 
possibly more sophisticated way around this problem is 
to employ a screened exchange instead of the full one£i22 

The outcome of the electronic structure methods intro- 
duced above is summarized in Tab. [J where we compare 
two electronic phases of MnO in the rocksalt structure 
at equilibrium lattice constant: the nonmagnetic (NM) 
state and the state with antiferromagnetic (AFM) order- 
ing of adjacent Mn (lll)-planes, i.e., the so-called type-II 
antiferromagnet. The AFM phase represents the ambient 
pressure ground state. In the case of the LAPW method, 
the basis size and the density of fc-point mesh was con- 
verged down to 1 or 2 mRy (~ 0.03 eV) in obtained total 
energies per formula unit. The Gaussian basis set for the 
LCAO methods was chosen large enough to reproduce 
the (nonrelativistic) LAPW/GGA results (see Tab-HJ). 

Since we utilized only small-core pseudopotentials 
we expected only minor discrepancies between the all- 
electron and PP calculations. However, this assumption 
turned out not to be entirely correct. The results listed 
in Tab. U show that both explicit presence of the core 
and the method which is employed do matter. The differ- 
ence between the all-electron and pseudopotential tech- 
niques varies from +0.6 eV/MnO within the PBE-GGA 
functional to —0.6 eV/MnO in the Hartree-Fock method 
combined with the STU PP. While contributions of this 
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magnitude might not be always relevant, for the purpose 
of comparing energies of different magnetic phases or in 
high pressure transitions such a difference turns out to 
be very significant. 

In addition to the antiferromagnetic and nonmagnetic 
Bl phases we included into our investigation also the fer- 
romagnetically (FM) ordered Bf structure. This time the 
energy difference A = Eafm — Efm stays virtually the 
same irrespective of the the core treatment used. This 
indicates that the sensitivity of A to the description of 
the core electrons is related to difference in magnetic mo- 
ments of the two involved states — the local moments in 
FM and AFM phases are almost identical, approximately 
4.5 fis (see Tab. HJ, whereas they are by definition zero 
in the NM state. More detailed evidence pointing in this 
direction is presented within the high pressure study bel- 
low. 

Several comparisons of all-electron and small-core 
pseudopotential data can be found in the literature in 
connection with the so-called non-linear core correction 
(NLCC)J^ LDA studies of 3d atoms 2 ^ and transition 
metal bearing diatomic molecules 24 show that without 
the NLCC the errors in energetics, associated with small- 
core pseudopotentials, are generally of the order of tenths 
of eV, which is in agreement with our observations made 
in solids. The non- linear core correction is demon- 
strated to be capable of removing a large portion of these 
errors. 23 - 24 
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TABLE II: Dissociation energy D e of MnO molecule at a bond 
length 1.64 A calculated with restricted open shell (ROHF) 
and unrestricted (UHF) Hartree-Fock methods. Both the 
molecule and Mn atom are in a high-spin configuration with 
spin multiplicity 6. 
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TABLE III: Energy difference A = Ei s — Ehs between high- 
spin (cr 1 5 1 5 1 w 1 jr 1 , multiplicity 6) and low-spin (ct 1 ^ 2 ^ 2 , mul- 
tiplicity 2) configurations of MnO molecule calculated with 
Hartree-Fock methods. Bond length is 1.64 A as in Tab. HTl 



We have decided not the include the NLCC in our 
study for the following reason. The NLCC concept, cor- 
recting for linearization of core-valence correlations in- 
herent to the standard norm-conserving pseudopotential 
technique, is not systematically transferable to methods 
using exact exchange, i.e., to hybrid exchange-correlation 
functionals, or to many-body methods such as quantum 
Monte Carlo. For these approaches it is therefore im- 
portant to assess how accurately the actual pseudopo- 
tential Hamiltonian represents the physical ion. Indeed, 
our study is a test whether this appears to be true for a 
variety of magnetic states, range of pressures and differ- 
ent methods. In the next section we provide additional 
benchmarks for MnO molecule and Mn atom, which 
clearly delineate the consistency of the used PPs and also 
suggest what is the reason for the discrepancies we found. 
Although the discussion of pseudopotential constructions 
is out of the scope of this paper we just note that the core- 
valence partitioning is more straightforward within the 
Dirac-Fock approximation than within the DFT due to 
the rather complicated nonlinear dependence of the DFT 
exchange-correlation functional on the density. In fact, 
one set of the Dirac-Fock pseudopotentials we employed 
has been used in high-accuracy correlated calculations by 
quantum-chemical methods in a routine manner without 
additional corrections and with satisfactory results for 
equilibrium properties of molecular systemsi 10 ' 11 ! 25 



III. TESTING THE PSEUDOPOTENTIALS IN 
Mn ATOM AND MnO MOLECULE 

The results presented in the preceding paragraph 
might appear somewhat surprising, since very small er- 
rors generated by small-core PPs have been observed in 
isolated atoms j 10 i 18 In particular, our Hartree-Fock tests 
of atomic s — ► d transfer energies, i.e., differences be- 
tween high-spin and low-spin states in the isolated atom 
(not explicitly shown), turned out the same no mat- 
ter whether all-electron or pseudopotential Hamiltonian 
has been employed. In addition, our previous calcula- 
tions within high-accuracy quantum Monte Carlo repro- 
duced the experimental s — > d transfer energy within 0.1 
eVi 1 ^ This suggests that the norm-conserving Dirac-Fock 
construction provides accurate PP Hamiltonian, which 
represents the physical ion within both the mean-field 
and correlated many-body methods with accuracy better 
than 0.1 eV. A question then arises where is the origin of 
the discrepancies reported in the previous section. Is it 
something related solely to the solid state environment or 
can it be detected as soon as chemical bonds are formed? 
The LDA s tudy of molecules containing transition metal 
atoms, Ref. [2J, suggests that the second statement is in- 
deed true. To verify this we performed several Hartree- 
Fock calculations to demonstrate that such a conclusion 
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is not only LDA specific. 

Tables [TT] and IIIII show results obtained using the 
GAMESS code-- with large uncontracted basis sets for 
pseudoatoms and with natural atomic orbital o 27 ^ 28 for 
all-electron atoms. Whenever only same spin states are 
involved, as is the case of the dissociation energy D e 
(Tab. |TT| of the ground state of MnO molecule, the all- 
electron and pseudopotential results agree quite well. Es- 
pecially when YN pseudopotential is used, the difference 
jjaii _ jjPP f a n s q U ite clearly bellow resolution of the ba- 
sis sets employed. In addition, tests done with highly cor- 
related multi-determinant Slater- Jastrow trial wave func- 
tions and with fixed-node diffusion Monte Carlo repro- 
duced the experimental binding energy within 0.1 eV*ii 
Therefore, even cohesive properties of the PP Hamilto- 
nian look correct both within mean-field and many-body 
correlated methods. 

A small but visible inconsistency between all-electron 
and pseudopotential approaches suddenly appears if 
one considers a high-spin to low-spin transition in 
the molecule (Tab. IIIip . Electronic configurations are 

S^utt^tt^ for the high-spin and a 1 ^ 2 ^ 2 for the low- 
spin state. Comparison of quantities A a " — A pp be- 
tween Tabs. HI and IIIII shows that the discrepancies in the 
solid are almost three times larger than in the molecule. 
This is understandable since the Mn bonding in solid is 
enhanced by the octahedral arrangement of the oxygen 
atoms. Trying to estimate how much stronger this effect 
would be from the size of cohesions one finds that the 
ratio of (experimental) binding energie o 29 ' 30 of the two 
systems is 9.5eV/3.8eV = 2.5 . The key question which 
therefore remains to be answered is: where is this effect 
coming from? We believe that the bonding environment 
is the key consideration for resolving this issue. In an iso- 
lated atom the d electrons have freedom to expand away 
from the nucleus and indeed one finds that the atomic 
d orbitals for low-spin states are very diffuse. In the 
presence of chemical bond, say, in a molecule, such an 
expansion is not favorable, since bonding tends to local- 
ize the states even for doubly occupied d channels. In a 
solid this is even more pronounced as the Mn atom is en- 
capsulated in the octahedron of oxygens. This enhances 
the amount of charge in the d states inside the core. The 
cause for the spin related discrepancies then becomes ap- 
parent — it is the bonding (or chemical) "pressure" . 



IV. SOLID AT HIGH PRESSURES (GGA) 

In order to examine the high pressure behavior of 
MnO, namely the Bl — > B8 transition, we have carried 
out a series of LAPW/GGA calculations^ for the type- 
II antiferromagnetic phase in the Bl structure and for 
the B8 structure with ferromagnetic as well as antiferro- 
magnetic ordering of hexagonal Mn planes.— The critical 
pressure P c and volume change at the Bl — > B8 transition 
are listed on the first two lines of Tab. IIVI The arrange- 
ment of local moments is found to change from AFM to 
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TABLE IV: Critical pressures for the Bl AFM B8 FM 
phase transition as calculated with various methods. Also in- 
dicated are volumes (per formula unit) and the local magnetic 
moments on Mn atoms for both phases at the transition. 



FM across the transition. The choice of the magnetic 
ordering in the high pressure phase is not well resolved, 
however, since the energy difference between B8 AFM 
and B8 FM phases at the relevant volumes is compara- 
ble to our 2 mRy precision estimate. The Fermi level lies 
in the gap between bands of the Kohn-Sham eigenvalues 
from ambient pressure all the way up to the transition 
in the case of the Bl AFM phase, whereas it is located 
inside a band for both studied B8 phases above and also 
well bellow P c . These qualitative features do not change 
if full cores or pseudopotentials are used. The variation in 
quantitative aspects is, on the other hand, quite sizeable. 
For instance, the difference in the transition pressure P c 
is as large as 25 GPa. 

The GGA predictions are rather far from experimen- 
tal P c around 100 GPaZ*^ and from a non-magnetic high 
pressure phase 8 observed at room temperature (our cal- 
culations would support a paramagnetic phase at this 
temperature, at best). This is not surprising as the GGA 
is not expected to approximate the electronic structure 
of MnO very accurately. At the same time, the pressure 
P c = 67 GPa resulting from our PP calculations is in 
a good agreement with results of Fang et al£^ who ob- 
tained P c w 60 GPa with a plane-wave pseudopotential 
approach. The remaining few GPa might originate in the 
full structure optimization performed in Ref. l3ll . For the 
sake of simplicity, we have neglected the rhombohedral 
distortion in the Bl AFM phase and have not optimized 
the c/a ratio in the B8 phases. Instead we used a con- 
stant value c/a — 2.0 consistent with outcome of Ref. l3ll 
and slightly smaller than c/a sa 2.1 suggested by x-ray 
spectroscopy.— 

Figures Q] and [2] relate the discrepancy between the 
two distinct core electron treatments to magnetic mo- 
ments formed on manganese atoms. The quantity we use 
to measure the difference between all-electron and pseu- 
dopotential total energies is defined as 

e(y) = E all (V)-E pp (V) - [E all (V ) -E PP {V Q )] , (2) 

where Vb stands for a reference volume. We use C = 
21.7 A 3 per formula unit in the Bl phase (Fig. Q} and 
Vb = 23.4 A 3 per formula unit in the B8 phase (Fig. [2]). 
In the rocksalt structure, the difference e remains almost 
constant as long as the spin polarization does not change, 
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FIG. 1: (color online) Difference between all-electron and 
pseudopotential total energies, Eq. ([2]), as calculated for the 
Bl AFM structure within LAPW/GGA method (bullets). 
Also shown is magnetic moment on the Mn atom with all- 
electron core (full squares) and YN pseudopotential (empty 
squares). The slight enhancement of local moments when 
pseudopotentials are introduced, visible also in Tabs.Hland llVI 
is consistent with other reports . 13 ' 32 Volumes at the transition 
to the B8 FM phase are indicated as well. 
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FIG. 3: (color online) Magnetic moment on Mn atom in Bl 
AFM (empty squares) and B8 FM (full squares) phases as 
calculated using the PBEOio functional. Values are taken 
from YN PP calculation, corresponding all-electron data are 
virtually indistinguishable. Transition volumes are shown as 
well. 



together. 




FIG. 2: (color online) Same quantities as in Fig. [T]but for B8 
FM structure. 

and jumps up at the high-spin to low-spin transition. 6 
The main reason for the 25 GPa inconsistency in the 
critical pressure P c lies in the B8 phase, however, since 
rather high pressure is necessary for the high-spin to low- 
spin transition in the rocksalt structure. In the B8 phase 
the magnetic moment decreases gradually with volume 
compression and this reduction is accompanied with a 
gradual change of the difference e that reaches an ap- 
preciable value at the point of the phase transition (see 
Fig. [2]) • The fact that the variation of the total energy 
difference e apparently follows the evolution of spin polar- 
ization suggests that some type of spin-dependent pseu- 
dopotentials, possibly those introduced in Ref. 33, might 
bring the all-electron and pseudopotential results closer 



V. SOLID AT HIGH PRESSURES (HYBRID 
FUNCTIONAL) 

To evaluate how the exact exchange influences the crit- 
ical pressure of the Bl — > B8 transition we repeated the 
calculations with the PBEOio functional as shown on the 
last two lines of Tablc HVl We can see that already a small 
admixture of exact exchange moves the critical pressure 
P c way above the GGA prediction, just a little higher 
than experimentally observed. Corresponding volumes 
of the phases above and bellow the transition are also 
very close to the experimental results^ We would like 
to emphasize, however, that this study is targeted on 
methodology rather than on precise comparison with ex- 
periments since the latter goal would certainly require 
detailed geometry optimizations. The overall picture of 
the transition within PBEOio stays the same as with the 
GGA, the rocksalt phase is an insulator, the B8 phase is 
a metal. The magnetic moments collapse from 4.4 fis to 
2.0 /is as the MnO is compressed across the transition. 
(Note that there are slight inconsistencies in comparisons 
of LAPW and LCAO local magnetic moments due to 
their different definitions.) The evolution of the mag- 
netic moments with decreasing volume is shown in Fig. [3] 
The tendency of the Fock exchange towards electron lo- 
calization is clearly visible, the decrease in moments is 
delayed to a lower volume in both Bl and B8 phases as 
compared to the GGA results displayed in Figs. [1] and [5] 
In the rocksalt structure we actually see no high-spin to 
low-spin transition in the interval of examined volumes. 

The discrepancy between all-electron and pseudopo- 
tential approaches is greatly reduced when a hybrid func- 
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tional is used. Following the discussion of the GGA re- 
sults one might expect that this reduction is accompanied 
by limited spin changes. However, the Fig. [3] demon- 
strates that this not the case and therefore we conclude 
that the core treated by PBEOio is more inert and is 
less sensitive to changes of the local moment formed by 
valence electrons. 



VI. CONCLUSIONS 

Our findings complement several other reports found 
in the literature indicating that application of pseudopo- 
tential based techniques to transition metals and their 
compounds is a delicate task. 14 In particular, we can 
mention discrepancies in magnetic properties of the vana- 
dium surfac o 34 i 35 ' 36 or inconsistencies in predictions of 
the ground state of Ce20 j 37138 found between pseudopo- 
tential and all-electron methods. In these studies, how- 
ever, the large-core (Ar) PPs were utilized. The diffi- 
culties with pseudopotentials are then attributed to the 
semi-core 3s and 3p states that have too large overlap 
with the valence 4s and 3d electrons to be accurately 
represented by a "rigid core". In the same time, the 
semi-core states are usually too localized to be efficiently 
included in the valence spaced 3 - This type of argument, 
however, cannot explain the variations in the Bl — > B8 
transition pressure presented in this paper, since we treat 
the 3s and 3p states as valence degrees of freedom. Our 
observations suggest that one should carefully examine 
transferability of pseudopotentials between different spin 
states even in the case of small-core PPs, especially if 
the phenomena under investigation involve variations in 
local magnetic moments. 

Overall implications of our results are rather straight- 
forward to explain. The 3d states have no counterparts of 
the same symmetry in the core and therefore an apprecia- 
ble amount of the d charge density is localized in the core 
region. Consequently, changes in the local moments can 
influence the core states and vice versa. As we have ver- 
ified, this is not detectable in isolated atoms, where the 



all-electron and pseudopotential results agree very well, 
nor in equilibrium ground state cohesion energies, which 
do not involve dramatic changes in the local moment on 
the Mn atom. This is true not only for the mean-field 
methods employed in this study, but also for highly ac- 
curate quantum Monte Carlo tests done previously ! 10 ! 11 
These results suggest that the PP Hamiltonian is con- 
sistent and accurate. The discrepancy appears and be- 
comes visible when the high- spin /low- spin differences be- 
come involved in bonds. The localization of the d states 
due to the chemical environment pushes larger fraction 
of the charge into the core and inaccuracies of various 
methods become more pronounced. The direction of the 
biases, which we found, is also not very difficult to un- 
derstand. The LDA and GGA functionals are well-known 
to incline towards more smooth and homogeneous densi- 
ties with exaggerated softening of electron-electron repul- 
sion. On the opposite side is the Hartree-Fock method, 
which tends to localize the states since that enhances the 
exchange contributions (the only mechanism, through 
which HF can decrease the energy). The key point is 
that these biases affect all degrees of freedom of the given 
system, including, in all-electron calculations, also the 
core states. Therefore, the cores are "softened" in the 
LDA/GGA approaches and "tightened" in the HF ap- 
proximation. In the PP case the core is rigid and there- 
fore the bias from core states is absent that leads to seem- 
ingly better (closer to experiment) results for the transi- 
tion pressures. Finally, these results could be of interest 
not only for high-pressure solids but also for other tran- 
sition metal systems with significant magnetic moment 
changes such as biomolecular reaction centers or cataly- 
sis. 
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